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Abstract 

The Smodels system implements the stable model se- 
mantics for normal logic programs. It handles a sub- 
class of programs which contain no function symbols 
and are domain-restricted but supports extensions in- 
cluding built-in functions as well as cardinality and 
weight constraints. On top of this core engine more 
involved systems can be built. As an example, we have 
implemented total and partial stable model computa- 
tion for disjunctive logic programs. An interesting ap- 
plication method is based on answer set programming, 
i.e., encoding an application problem as a set of rules 
so that its solutions are captured by the stable models 
of the rules. Smodels has been applied to a number 
of areas including planning, model checking, reachabil- 
ity analysis, product configuration, dynamic constraint 
satisfaction, and feature interaction. 



General Information 

The Smodels system is written in C++ and the source 
code, test cases and documentation are available at 
http://www.tcs.hut.fi/Software/smodels/. In or- 
der to compile the system a C++ compiler is needed 
as well as other standard tools such as make, tar, bi- 
son, and perl. The system has been developed under 
Linux and should work as is on any platform having the 
appropriate GNU tools installed. It has been used on 
a wide range of hardware (PC/SPARC/ Alpha) mostly 
running Unix. The total number of lines of codes is 
about 20000. 

Description of the System 

The Smodels system implements the stable model 
semantics for normal logic programs extended by 
built-in functions as well as cardinality and weight 
constraints for domain-restricted programs. In 
this section we briefly discuss the syntax, im- 
plementation techniques and use of the system. 
More information can be found at the home page 
http : //www. tcs .hut . f i/Sof tware/smodels/. 



*The work has been funded by the Academy of Finland 
(Project 43963) and by Helsinki Graduate School in Com- 
puter Science and Engineering. 



As input the Smodels system takes logic program 
rules basically in Prolog style syntax: 

ancestor (X,Y) :- ancestor(X,Z) , 

parent (Z, Y) , 

person(X) . 
ancestor (X,Y) :- parent(X,Y). 
son(X,Y) :- parent(Y,X), male(X). 
daughter(X.Y) :- parent (Y,X) , female(X). 
person(X) :- male(X). 
person(X) :- female(X). 
parent (jack, jill). parent (joan, jack). 
male(jack). f emale ( j ill) . female (joan) . 

However, in order to support efficient implementation 
techniques and extensions the programs are required 
to be domain-restricted where the idea is the follow- 
ing. No proper function symbols are allowed (but we 
do allow built-in functions) and the predicate symbols 
in the program are divided into two classes, domain 
predicates and non-domain predicates. Domain pred- 
icates are predicates that are defined non-recursively. 
In the program above all predicates except ancestor 
are domain predicates. The predicate ancestor is not 
a domain predicate because it depends recursively on 
itself. 

The main intuition of domain predicates is that they 
are used to define the set of terms over which the vari- 
ables range in each rule of a program P. All rules of 
P have to be domain-restricted in the sense that every 
variable in a rule must appear in a domain predicate 
which appears positively in the rule body. For instance, 
in the first rule of the program above all variables ap- 
pear in domain predicates parent and person in the 
body of the rule. 

In addition to normal logic program rules, Smodels 
supports rules with cardinality and weight constraints. 
The idea is that, e.g., a cardinality constraint 

1 { a.b.not c } 2 

holds in a stable model if at the least 1 but at most 2 of 
the literals in the constraint are satisfied in the model 
and a weight constraint 

10 [ a=10, b=10, not c=10 ] 20 



holds if the sum of weights of the literals satisfied in the 
model is between 10 and 20 (inclusive). With built-in 
functions for integer arithmetic (included in the sys- 
tem), these kinds of rules allow compact and fairly 
straightforward encodings of many interesting prob- 
lems. For example, the N queens problem can be cap- 
tured using rules as a program queens . lp as follows: 

1 { q(X,Y):d(X) } 1 :- d(Y) . 
1 { q(X,Y):d(Y) } 1 :- d(X) . 
:- d(X), d(Y), d(Xl), d(Yl), 

q(X,Y), q(Xl,Yl), 

X != XI, Y != Yl, 

abs(X - XI) == abs(Y - Yl) . 
d(l. .n) . 

where d ( 1 . . n) is a domain predicate giving the dimen- 
sion of the board (from 1 to integer n which can be 
specified during run time). The first rule says that for 
each row y, a stable model contains exactly one atom 
q{x, y) where for x, d(x) holds and similarly in the sec- 
ond rule for all columns. The third rule is an integrity 
constraint saying that there cannot be two queens on 
the same diagonal. Now each stable model corresponds 
to a legal configuration of n queens onanxn board, i.e., 
q(x,y) is in a stable model iff (x,y) is a legal position 
for a queen. 

Stable models of a domain-restricted logic program 
with variables are computed in three stages. First, the 
program is transformed into a ground program with- 
out variables. Second, the rules of the ground program 
are translated into primitive rules, and third, a stable 
model is computed using a Davis-Putnam like proce- 
dure ( Simons 1999 ). The first two stages have been 
implemented in the program lparse, which functions 
as a front end to smodels which in turn implements 
the third stage. 

In the first stage lparse automatically determines 
the domain predicates and then using database tech- 
niques evaluates the domain predicates and creates a 
ground program which has exactly the same stable 
models as the original program with variabl es. Then 



the rules are compiled in to primitive rules ( Niemela 
Sim 



ms, fc Soinincn 1999 ) 



The smodels procedure is a Davis-Putnam like back- 
tracking search procedure that finds the stable models 
of a set of primitive rules by assigning truth values to 
the atoms of the program. Moreover, it uses the prop- 
erties of the stable model semantics to infer and prop- 
agate additional truth values. Since the procedure is 
in effect traversing a binary search tree, the number of 
nodes in the search space is in the worst case on the 
order of 2", where n can be taken to be the number of 
atoms that appear in a constraint in a head of a rule or 
that appear as a negative literal in a recursive loop of 
the program. 

Hence, in order to compute stable models, one uses 
the two programs lparse, which translates logic pro- 
grams into an internal format, and smodels, which 
computes the models, see Figure [l]. 



For instance, a solution to the 8 queens problem given 
the program queens . lp above, would be typically com- 
puted by a command line: 

lparse -c n=8 -d none queens. lp I smodels 

where -c n=8 option instructs to use the value 8 for the 
constant n and -d none option instructs to remove all 
domain predicates from the rules as soon as they have 
been evaluated. The command line produces output: 

Answer: 1 

Stable Model: q(4,l) q(2,2) q(7,3) q(5,4) 
q(l,5) q(8,6) q(6,7) q(3,8) 

A Further Example 

The graph coloring problem may be encoded to a 
program ncolor . lp using the following two Smodels 
rules: 

1 { col(N, C) : color(C) } 1 :- node(N). 

:- col(X, C), col(Y,C), edge(X.Y), color(C) . 

Here the predicate col(N,C) denotes that the color of 
the node N is C. The first rule states that each node has 
exactly one color and the second rule states that two 
adjacent nodes may not have the same color. Suppose 
we have a fully connected three node graph and we want 
to find all 3-colorings of it where the first node is colored 
red. We can encode the problem instance to a program 
graph. lp: 

node(a ; b; c) . 

edge(a,b). edge(a,c). edge(b,c). 
color(red ; green ; blue) . 
compute { col (a, red) }. 

The first two lines define the graph and the third line 
defines the three colors. The compute statement tells 
Smodels that we are interested only in those models 
where col (a, red) is true, i.e., the node a has color 
red. We can find all stable models that satisfy the 
compute statement with the command line: 

lparse -d none ncolor. lp graph. lp I smodels 

This input produces the output: 
Answer : 1 

Stable Model: col(c,blue) col(b, green) 

col(a,red) 

Answer : 2 

Stable Model: col(c, green) col(b,blue) 
col(a,red) 

More information about the syntax and use of the 
system can be found in the lparse user's manual at 
http : //www. tcs .hut . f i/Sof tware/sm odels/lparse / 
and about the implementation in ( Bimons 1999] ; 
Nicmcla, Simons, fc Soinincn 1999[ ). 



Applying the System 
Methodology 

An interesting application methodology f or Smod- 
ELS is based on answer set programming (Marek & 
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Figure 1: Overall Architecture 



Truszczynski 1999; Niemela 1999 ) which has emerged 
as a viable approach to declarative logic-based knowl- 
edge representation. It is based on the stable model 
semantics of logic programs and can be seen as a novel 
form of constraint programming where constraints are 
expressed as rules. The underlying idea is to encode an 
application problem using logic program rules so that 
the solutions to the problem are captured by the sta- 
ble models of the rules. The solution of the N queens 



prol plcm in the previous section illustrates nicely main 



Users and Usability 

The basic semantics of Smodels rules seems to be sim- 
ple enough that it can be explained to a non-expert or a 
student in a relatively short amount of time sufficiently 
so that one can start using Smodels. 

Smodels has b een employed in a number of areas in- 
cluding planni n g flDimopoulos, Ncbcl, fc Kochlcr 19971: 
Niemela 1999]; Lifschitz 1999|), model checking (Liu 



ideas of answer set programming. 

Sp ecifics 

It 



Ramakrishnan, fe Smolka 1998), rea chability analy 
sis ( |Hcljanko 1999a; [Heljanko 1999b), product con - 
figuration ( goininen fc Niemela 1999]; [Byrjanen 1999 \ 



dynamic co nstraint satisfaction 
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Gelle 1995;), feature interaction QAccorsi, Areces, fc dc 



Soinincn, Niemela^ & 



Rijke 1999), and logica l cryptanalysis (Hictalahti, Mas 



Thia makes it much easier to develop applications be- 



sacci, fc Niemela 2000) 



cause one does not have to worry too much about 
the internal implementation-specific aspects of the sys- 
tem. Hence, the system is relatively easy to learn to 
use. On the other hand, declarative semantics provides 
much more flexibility in developing implementation ap- 
proaches and in optimizing different parts of the imple- 
mentation. We have taken advantage of this and de- 
veloped methods which are substantially different from 
usual implementation methods for logic programming 
(Prolog) but still work efficiently in new kinds of appli- 
cations where Prolog style systems are not appropriate. 

Smodels implements the stable model semantics for 
range-restricted function free normal programs. It can 
also compute well-founded models for these programs. 
Smodels supports built-in functions, e.g., for integer 
arithmetic. Basic Smodels extends normal lo gic pro- 
grams with cardinality and weight constra ints ( Simous| 
199{fl; Niemela, Simons, & Soinincn 1999). On top of 



this core engine more involved systems can be built. 
As an example, we have implemented total and par- 
tial sta ble model computation for disjunctive logic pro- 
grams (Janhuncn et al. 200C). 

The semantics for logic programs with cardinality 
and weight constraints supported by the core engine 
is an interesting compromise: it is rather simple to 
learn, its complexity stays in NP in the ground case 
(like for propositional logic) but it seems strictly more 
exp ressive than propositional logic (or other standard 



cons traint satisfaction formalisms). 



With our work on Smodels wc hope to demonstrate 
that nonmonotonic reasoning techniques are useful con- 
ceptual tools as well as bring computational advantages 
which can lead to new interesting applications and to 
developments of novel implementation techniques and 
tools. 



In order to make Smodels more flexible an API has 
been added to Iparse and smodels. Hence, they can 
be used through the API and embedded into a C/C++ 
program as libraries. Furthermore, it is possible to de- 
fine new built-in functions for the front-end Iparse in 
order to accommodate new applications. These new 
functions are written in C/C++ and they are dynami- 
cally linked to Iparse when needed. 

As a further usability feature, Iparse performs sim- 
ple analysis of the program and warns about constructs 
that are often erroneous. For example, Iparse detects 
if a variable is accidentally mistyped as a constant or 
vice versa. These warnings can be enabled with com- 
mand line options. 

Evaluating the System 
Benchmarks 

We have compiled quite a large collection of families of 
benchmarks that we use to evaluate new developments 
and improvements in our system and which we can 
employ to compare our system to other compet- 
ing ones and which can be used also by others (see, 
http : //www . tcs .hut . f i/Sof tware/smodels/tests/). 
Similar testing methodology (e.g., generating test cases 
from graph problems) has also been used b y the groups 
in Kentucky (DcReS/TheoryBase system ( ICholcwiiiski 



et al. 1999)) and in Vienna (dlv system (Eitcr et al 



1998)). These kinds of tests can be used for measuring 



the base level performance and to compare different 
systems. However, it is unclear what is a portable way 
of representing such benchmarks. We use Prolog style 
syntax which is quite generally accepted but it seems 
that each system has its specialties. 



Comparison 

Smodels can already compete with special purpose 
systems and we have even cases where it outperforms 
commercial top edge tools, e.g., in a verification ap- 
plication it has performance better than that of one 
of the most effic ient commercial integer programming 



tools (CPLEX) (Hcljanko 1999b) 



Problem Size 

We believe that Smodels is no longer a mere prototype 
but it can handle realistic size problems. We have ap- 
plications where programs with hundreds of thousands 
of n on-trivial ground rules are trea ted efficiently, see 
e.g. flNiemela 1999j jHeljanko 1999b|) . 
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